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ABSTRACT 
The deep basins of the world oceans are filled with waters of 
very low temperature. A fundamental problem of oceanography concerns 
the nature of the sources of that water: their locations and temporal 
variations of strength. A stability analysis of temperature and salinity 
data recently taken in the Norwegian sea shows that bottom water probably 


does not form in that area. PN ae 
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dy, Introduction 

The deep basins of most of the world oceans are filled with waters 
of very low temperature. The source of these cold bottom waters has 
been the subject of considerable research in oceanography. Nansen, 
in 1912, formulated what has become the accepted theory for the forma- 
tion of these waters. His theory, stated briefly, is as follows. 
During the cold season, the salinity of the surface layer in the high 
latitudes is increased either through advection of more saline water 
from lower latitudes or through freezing of the surface water in situ. 
Concurrent with or subsequent to this salinity increase the temperature 
of the surface water is lowered by heat loss to the atmosphere, through 
the processes of radiation, convection and conduction. The combina- 
tion of increased salinity and decreased temperature produces surface 
water which is more dense than the stratum of water immediately below 
the surface layer. This more dense water will sink, mixing with the 
underlying waters and producing a thin layer of homogeneous water. As 
the process of formation proceeds the thickness of this homogeneous 
layer will increase, and if enough dense surface water forms the layer 
will ultimately extend from the surface to the bottom. If the process 
of dense water formation is not sufficiently intense or if it does not 
continue for a sufficient length of time the homogeneous layer will not 
extend to the bottom but will terminate at an intermediate depth and be 
destroyed by lateral mixing at the conclusion of the surface formation 
period for the dense water. 

Since the formation process for bottom water depends upon inter- 


action between the sea and the atmosphere it is not likely that the 


process proceeds at a uniform rate, but rather that it occurs at greater 
or lesser rates depending on variations in the upward transfer of heat 
and in the surface salinity. Variations in wind, precipitation, air 
temperature, etc., should vary the rate of bottom water production over 
short intervals of time, perhaps of the order of days or even hours. The 
result would not be a column of water having a uniform density from sur- 
face to bottom, but rather a column having vertical density variations 
corresponding to periods of formation of waters of different density at 
the surface. During a relatively long period of formation sufficient 
water of high density might be produced so that a stratum could retain 
its characteristics long enough to enable it to reach great depths with- 
out complete loss of identity through mixing, and thus be detectable 

as inhomogeneity in the density distribution. 

Static vertical stability in a fluid is a measure of the restoring 
force on a fluid element when it is given a small adiabatic vertical 
displacement from its undisturbed position. This force is proportional 
to the difference between the density of the displaced element and the 
density of its surroundings. In a stably stratified fluid an element 
displaced downward will be less dense than its surroundings and the 
resulting buoyant force will tend to return the element to its former 
location. In a neutrally stratified fluid the displaced element will 
have the same density as its surroundings after displacement and will 
experience no buoyant force. In an unstably stratified fluid the ele- 
ment, on being displaced downward, will have a density greater than that 
of its surroundings and will experience a negative buoyant force tend- 
ing to displace the element farther from its former location. 
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In an area of bottom water formation where atmospheric conditions 
vary with time, the stability of the surface layer should vary from 
stable through neutral to unstable, depending upon the heat and salt 
content. After the beginning of formation of actual bottom water (after 
the initial formation of the surface-to-bottom homogeneous column) the 
stability of the deeper columns should be essentially neutral. Finite 
strata of unstable water might be interspersed from time to time, if 
there have been periods of extensive formation of very dense water. 

According to Sverdrup et al. (1942), the favored areas for forma- 
tion of bottom water in the North Atlantic Ocean are the Labrador, 
Irminger, and Norwegian Seas. An extensive time series of hydrographic 
station data is available for Weather Station Ship Mike (Mosby undated). 
Ship Mike is located in the southeastern Norwegian Sea at latitude 66° 
N. longitude 02° E. Figure (1) shows the bathymetry of the southern 
Norwegian Sea and the location of Ship Mike. These data were analyzed, 
through stability calculations, to determine whether or not bottom 


water formed at that location and, if so, the temporal variation of 


the process. 


2. The Data 

The data used for this investigation consist of 1253 hydrographic 
stations taken at Ship Mike during the period January 1949 through 
December 1953 (Mosby undated). It was obtained in machine form (punched 
cards) directly from the Geophysical Institute, University of Bergen, 
through the cooperation of Professor Hakon Mosby. Observations were 
made using four different combinations of depths for the Nansen bottles 
and reversing thermometers. The four combinations and the number of 


observations made using each are shown in Table l. 


TABLE 1 


SAMPLE DEPTHS FOR THE FOUR TYPES OF OBSERVATIONS TAKEN AT SHIP MIKE AND 
THE NUMBER OF EACH TYPE OF OBSERVATION TAKEN 


1949-1953 
Type : 
Depths aes ee mo ss 
(meters) 0 0 0 0 
10 10 10 50 
25 25 25 150 
50 50 50 300 
75 75 75 600 
100 100 100 1000 
150 150 150 
200 200 
300 400 
400 600 
500 1000 
600 
800 
1000 
1200 
1500 
2000 
235 391 231 396 Total number of each 


type 
The four combinations are fairly evenly distributed among the 


months. Occasionally a month passed with no type (a) observation but 


10 





for most months there is at least one and are as many as five observa- 
tions at depths greater than 1000 meters. 

Mosby (1959) has carefully evaluated the accuracies of the observa- 
tions. In general, standard oceanographic accuracies in temperature, 
salinity and depth were achieved. However, four to five percent of the 
salinity determinations contained errors greater than 0.01 parts per 
thousand. Observations having this reduced accuracy are not identified 
in the data. There are position errors in the location of the observa- 
tions, due presumably to the inability of the observing vessel to main- 
tain station in heavy weather. Most of the observations were taken 
within 15 to 20 nautical miles of the assigned station. With very few 
exceptions the observations were made in the morning hours, at about 


1000 local time. 


| 


3. The Stability Calculations 


Nansen's theory for the formation of bottom water has been pre- 
viously discussed. A primary requirement for such formation is the 
presence of a column of water, extending from the surface or near sur- 
face layer to the ocean floor, which exhibits small stability. That is, 
the column must have a nearly uniform density from top to bottom, ex- 
clusive of density changes due to adiabatic heating. The calculation 
of vertical stability was thus selected as the means for determining 
the occurrence of conditions suitable for the formation of bottom water 
at Station Mike. Vertical stability can be calculated using the formula 


of Hesselberg and Sverdrup (Sverdrup et al,1942): 

E = side + 20) = (Hh), (1) 

dS dh = Srl dh dh J, 

where c is density in Sg is salinity (parts per thousand), | is 
temperature, h is depth, and (ST/Ah ), represents the adiabatic 
temperature change with depth. Here, E is the net restoring force 
per unit volume and displacement, divided by the gravitational accelera- 
tion, 9 . This equation is inconvenient for digital computer calcula- 
tions. The terms dS/qh and d1TMh are evaluated directly from the 
data. The other terms are evaluated using extensive tabulated data. 
The inclusion of these tables lengthens the calculation time and reduces 
the amount of the memory core available for output data storage during 
the calculation process. The stability parameters developed by Pollack 
(1954) are much more convenient for computer use and were chosen for 
this investigation. 


Pollack's formulas are divided into three families. Calculations 


1 


made with the first family give stabilities in terms of force per unit 
mass per unit displacement, denoted by E.. Those made with the second 
family give stabilities in terms of potential energy per unit mass per 
unit displacement, ES Results using the third family give stabilities 
as a dimensionless ratio between the in situ density gradient and the 


adiabatic density gradient, The subscript (i) depends upon the 


ere 
method of measuring depth, which may be in terms of pressure (i = p), 

geopotential (i = G), or geometrical length (i = h ). Equations (2), 
(3), and (4) are examples of Pollack's formulas. 


21 8e- (29 ‘ 


et = at [ Se -(HF),| . 


ean = GE [CE “ 


Here, ( de/dh ds is the adiabatic density gradient. Since the 


speed of sound is given by c = \ (2 e/.4 ea , the above equations 


can be written in the more convenient forms: 


= HE 2 ° 


x ae Je _ et | 
2 ie ab ae (6) 


a 
= af 7 
Eeouny eS a ” 





13 


The dimensionless stability parameter, aye was selected for the 
stability calculations, and will hereafter be written as E- Values of 
E. greater than unity indicate stability, values equal to unity indicate 
neutral stability, and values less than unity indicate instability. 
Negative values of E indicate a decrease in density in situ with depth. 


Following Pollack, E is related to E by 


h 
a, 
and to E. by 
sf rt Ceicers (9) 
c 
A fairly accurate and simple approximation to equation (9) can be 
obtained. Assuming that ce. lenvene ; Pees 10° enVyeeer and 
ne = 2x Tere emaece equation (9) can be written as 
6 a para al 
E = = = a ee ae) (10) 


Stabilities calculated using Pollack's equation compare closely 
with those calculated using the Hesselberg-Sverdrup equation. Table 
2 illustrates stability calculations made using several of the equa- 
tions described above. The last two columns have been added to a similar 
table given by Pollack. The degree of approximation resulting from the 
use of equation (10) can be determined by comparison of the fourth and 


seventh columns of the table. 
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TABLE 2 


COMPARATIVE VALUES OF VARIOUS STABILITY PARAMETERS ON A SINGLE SET OF DATA 
(Note: we here deviate from standard oceanographic practice by using 
cgs units in E, EV» and E 


Depth Temp. Salinity 10°F 10’F, E 10°F from 
m a 8 Ppt Sec 2 : . eq: (10) 
« é 4 
Tee gm /cm 
0 19.20 36.87 
-440 -434 05 -475 
10 19.31 JO285 
-150 -167 07 -160 
25 19.34 30283 
- 13 - 15 v6 - 20 
50 19.24 30279 
610 604 ye 665 
75 18.65 36-79 
390 364 Le 425 
100 18.24 36.78 
38* 64 1.04 20 
150 17250 36.56 
270 242 Leo 280 
200 16.45 36.40 
160 160 1.30 150 
300 L432 36.02 
120 127 be 125 
400 13.08 iG Ba me a 
150 L152 1.34 170 
500 Die. 35.64 
130 134 ere) 145 
600 10.80 35.54 
100 107 Laz 115 
800 9709 goo 
89 96 | yt 105 
1000 8.01 i hoe 
84 86 1.20 100 
1200 reat | 35.42 
48 50 jE 60 
1400 6.40 a5255 
26% 29 L307 BS 
2000 ae. i os, 
11 14 Lae 15 
3000 2.84 eS 
8 Et 1.01 5 
4000 2.43 34.90 
i 3 99 -5 
5000 2240 34.90 


*Pollack recomputed these values using the formulas and tables of 
Hesselberg and Sverdrup 
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4. The Calculation of Variables 

Calculations: of E using equation (7) require values of both the 
speed of sound in sea water and the density of sea water in situ. 
These values are normally obtained from tabulated data (e.g. LaFond 
1951). Empirical formulas for the speed of sound (Kinsler and Frey 
1962) and for the density in situ (Eckart 1958) are available. These 
formulas eliminate the need for tabulated data and are thus more con- 
venient for use in computer calculations. 

The Kinsler-Frey formula for the speed of sound is 


3 
oe ee eee eS S I eee, Ss os 
(in) 
+ (1.39-0.01Q7T)(S-350) + 0.017 h 
where c is in m/sec, T in degrees centigrade, S in parts per thousand 
and h in meters. Results obtained using equation (11) are accurate to 


within one meter per second. 


Eckart's formula for the calculation of density in situ is 


era. £>) 
IO K = 19.995 t0.11gsT- 0.000 74ST (3625020, Coc. 
Poem SEAO + 38ST - 0. 275 T7274 23S 
No = 73020 (12) 
K = (e+ Po)(Ao—A). 


Here, P is measured in pnai/ che. T in degrees centigrade, S in parts per 
thousand, while © denotes absolute pressure in atmospheres. In the cal- 
culations the approximate expression p = 0.101 h was used. This approxi- 


mation introduces an error of about 5 atmospheres at a depth of 5000 meters 
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and may account for the slight instability in the 4000-5000 meter layer 
in Table 2. 

Eckart's formula, based upon almost all of the experimental deter- 
minations of the density of sea water, is valid for the ranges 0 to 40 
degrees centigrade temperature, 0 to 40 parts per thousand salinity, 
and 0 to 1000 atmospheres pressure. For these calculations the temper- 
Table 3 contains comparative 


ature range has been extended to -1.0°C. 


values of density calculated using LaFond's tables and Eckart's formula. 


TABLE 3 


DENSITY CALCULATIONS USING LAFOND'S TABLES AND ECKART'S FORMULA 





Depth Temp Salinity Tables Difference Eckart 
m = Ppt gm/em? gm /cm> gm/cm> 
200 6.35 35.19 1,02636 .00014 1.02852 
300 4.50 35.06 1.02918 . 00004 1.02914 
400 2.30 34.97 1.02983 . 00002 1.02981 
500 0-5 34.94 1.03039 00002 1.03041 
600 0.20 34.93 1.03092 .00004 1.03096 
800 -0.40 34592 LO3169 .00009 1.03198 

1000 -0.62 34.92 1.03284 .00012 1.03296 

1200 -0.75 34.92 1.03378 .00015 103293 

1500 -0.87 34.92 L703519 .00018 1203537 

2000 -0.95 34.92 L,03751 .00023 103774 


The data for 


Mosby (1959) 


at each level are apparent. 


depth, 


as 


"normal" 


1 
for Station Mike. 


temperature, and salinity are those described by 
The differences in density 


According to Eckart, there is no reason to 


ren determined his "normal" distribution of temperature and salinity 
with depth by finding the depth, on each of more than 200 hydrographic 
soundings which reached 2000 meters, of the integer isotherms and isohalines. 
These depths were then averaged and the resulting values plotted against 
the standard depths. In this way, according to Mosby, truer representa- 
tions of the average gradients of temperature and salinity were retained. 


Ly 


believe that the LaFond results are more accurate than those obtained 
using the Eckart formula. The significance of the differences and 
their relation to the calculation of stability will be discussed below. 

The stability of a column of water is calculated layer by layer. 
When evaluating equation(7) for EO the average speed of sound and the 
average density in each layer are used. There are two ways in which 
these average values can be determined. First, the values of tempera- 
ture, salinity, and depth for each level in a hydrographic sounding can 
be substituted into equations (11) and (12). The resulting values of 
the speed of sound and density can then be averaged and considered to 
represent the "average" sound speed and density within the layer. Second, 
the values of temperature, salinity, and depth can be averaged within 
each layer and the averages substituted into equations (11) and (12). 
These results also represent the "average'’ sound speed and density with- 
in the layer. Table 4 shows densities and speeds of sound calculated 
using the two approaches. The data are typical values of depth, tempera- 
ture, and salinity for Station Mike. The columns entitled Oy and ay 
represent the arithmetic averages of the values of e and C at the 
top and bottom of the layer. The columns headed O4¢h, av and Crehay 
represent the calculations made with equations (11) and (12), using the 
average values of T, S, and h within each layer. 

The question raised is one of linearity. Which are more linear 
with depth, the speed of sound and density or the temperature and 
salinity? An adequate answer to the question requires accurate continu- 
ous determinations of these parameters Simrleencouniyea For these 


: apt 
oon 


2 
The small discrepancies between the two methods are well within 
the maximum probable error due to other causes, as will be shown below. 
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calculations the values of temperature, salinity and depth were averaged, 
since computer calculations of simple arithmetic means are faster than 


calculations of quadratic and higher order terms. 


1 


te 2 ee 
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5. The Accuracy of E. 


The accuracy of E depends upon the accuracy of the calculation of 
the speed of sound, the accuracy of depth determination and the precision 
with which density in situ and its variation with depth can be obtained. 
The limit of accuracy of the sound speed calculation is one meter per 
second. Thermometric depths are, according to Von Arx (1962), accurate 
to within five meters at depths less than 1000 meters and to within 0.5% 
of the depth at greater depths. The effect of these errors on the ac- 
curacy of E can readily be examined. 

Assuming an absolute layer thickness error of §4) and an absolute 


sound speed error of dc , the calculated value of E. is 


(c+ ay do 





(Ee. ) (13) 
7 Pere dh+ fh 
where dh = hio-h, - Then the maximum error in E. is given by 
, Bc awa 
= -(€ = = —— >), 14 
nee) =) a =.} a a =} (14) 
Setting ¢h = 0.005 (h, +h,) and dc = 1 gives 
A E-(h,<) —- Ey, f = = 2.005 Chi th) : (15) 


eset) 3 ~ 


The maximum values for hy and h, at station Mike are 1500m and 2000m. 


A typical value for c is 1470 m/sec. Thus, 


KE. Chic) = CO4E,. (16) 
This result holds only for layers below 1000m , but the order of 
magnitude of AE. Chee) should be the same above this level. 


Zk 


The foregoing indicates that errors in Dh resulting from errors in 
the speed of sound calculation and thermometric depth determination will 
not exceed about 4% at Station Mike. 

Errors in the determination of density are less well known. The 
effect of these errors will now be examined. 

Disregarding the small systematic errors in the equation of state, 
the density errors are given by 
and (ep. = p + do, 
where Se represents the error in density determination. Then the 


error in E due to density determination is 


AE, Ce) = (E-).- ae = e-{ 2 dp + J { - (17) 
(doe CP). 


The committee on Oceanography, National Academy of Sciences (Vetter 
1959) stated that specific volume determinations(and hence, to a good 
approximation, density determinations) can be considered accurate to 
better than one in 10,000. Thus, taking e = 1 Eye de - 


0.001 ny eae and Je = 0.0001 ey ae 


PE eGe r= soraea ax (18) 


A 20% error in E. is possible in the stability calculations. 


Ze 


6. The Results of the Calculations 

The stability, EL was calculated for each layer of the 1253 hydro- 
graphic stations taken at Station Mike during the period January 1949 
through December 1953. The monthly means for E. were then calculated. 
The results are shown in figures 2, 3, 4, 5, and 6. 

Certain features are evident after an examination of these figures. 
The surface layer (0 to 25 meters) exhibits stability during the summer 
months and instability during the winter months except for the year 1953, 
and the duration of the surface instability is longer than the duration 
of the surface stability. The deeper waters (below about 800 meters) 
exhibit small stability with E. values only slightly greater than unity 
throughout the year. The striking feature, when considering the bottom 
water formation, is the presence and persistence of an intermediate 
layer of stable water. In the summer season this stable layer extends 
from 25 meters to about 600 meters. In the winter months it is somewhat 
more variable, but generally is located between 300 and 700 meters, and 
is never less than 200 meters thick. In the spring and fall the stable 
layer varies between the winter and summer distributions. 

The stable layer discussed above has been delineated using an ES 
value of 1.20 as the dividing line between stability and instability, 
in keeping with the probable accuracy of E That this is a reasonable 


choice can be further demonstrated. Recalling that E* is the potential 


h 
energy per unit mass and displacement, it follows that 
an . 2 
a dH = — 
ad (19) 


hy 
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where Vis the incident velocity a water element must have in order 
to pass downward through the layer limited by depths hy and h,- Since 
EX and E are related by 


* o* 
re een 


it follows that 
2 


\/* 7 3 CELy 1) Ch, - hia), (21) 


For the typical winter thickness of the stable layer, 200 meters, and 
speed of sound at the depth of this layer, 1468 m/sec, equation (21) 
yields for an E value of 1.20 an incident vertical velocity of about 
4.5 cm/sec. This represents the downward speed which a quantum of water 
would have to achieve to penetrate the stable layer and reach the bottom. 
Such large vertical velocities are unlikely to occur. 

The mean stability values have another striking feature. The mean 
depth of the stable layer is, in all cases, between 600 and 800 meters. 
Its exact location is unknown, due to the large spacing between observa- 
tions in the deep hydrographic casts. However, on six different occa- 
sions during the period 1951-1953 the depth of the base of this layer 
decreased sharply, moving upward to about 400 meters. This represents 
an increase in the thickness of the bottom-water layer of roughly 400 
meters. 

Also evident in the mean stability patterns are quanta of water 
having E values less than unity, indicating instability. These quanta 
appear from time to time during the winter and spring months, but in all 
cases they remain Hbece the stable layer discussed above. 
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The question arises concerning the stability distribution in the 
individual days of the period of observation. Figures 7 and 8 show the 
distribution of E. for a typical winter and summer month. The upper 
layers exhibit more variability than the monthly mean patterms, but the 
intermediate and deep layers closely resemble the monthly mean pattern. 
The individual deep observations show that the stable layer is never less 
than 200 meters thick. The value of the stable-layer E. in the individual 
observations falls below 1.20 only on six occasions and then to a value 
no smaller than 1.16. On each of these six occasions the thickness of 
the stable layer was at least 400 meters, which more than offsets the 
effect of the slight decrease in E. in the above velocity calculation. 

Examination of the individual observations during the periods when 
the base of the stable layer rises confirms the increase in thickness 
of the bottom waters. In all six cases there is a regular reduction in 
E. between 500 and 800 meters, from typical values of from 1.20 to 1.30 
to values 1.05 to 1.09. The reduction persists during the period of 
the thickness increase in the mean pattern and then returns to its normal 
value. 

Finally, the individual observations support the presence of un- 
stable water layers. However, due to the non-constant interval between 
observations, no conclusions can be drawn concerning the vertical move- 
ment of this water, except that it always appears to remain above the 
stable layer. 

The ability of the calculations to detect stability conditions 
associated with Nansen's theory for bottom water formation was tested 


in another manner, using data from an unpublished manuscript by Gladfelter 


Z 


(1964). These data were taken in the Greenland Sea, an area where 
bottom water is known to be formed. The vertical stability distribu- 


tions at two stations are shown in Table 5. 


TABLE 5 


STABILITY DISTRIBUTIONS IN THE GREENLAND SEA 


Station 1 Station 2 
Layer 74NOSW 7 3NOOW 

E E 

r r 

0-10 1.07 1.15 
10-20 0-05 8.21 
20-30 8.39 9.43 
30-50 3.80 3.06 
50-75 1.48 Te? 
75-100 ee el 1.34 
100-150 ye ee 0.99 
150-200 1.04 1.18 
200-250 1.04 1.02 
250-300 1.03 1.07 
300-400 1.06 1.03 
400-500 1.01 1.02 
500-600 1502 1.02 
600-800 1.03 1.02 
800-1000 1.01 1.04 
1000-1200 1.02 0.99 
1200-1500 120) 1.01 
1500-2000 1.01 1.01 
2000-2500 1.00 1.01 
2500-3000 1.01 1.00 


The results show clearly the stability characteristics for bottom 
water formation. With the exception of the surface layer both columns 
exhibit nearly neutral stability. Gladfelter's analysis of these sta- 
tions (and many more in the vicinity) indicates that this region is one 
of massive bottom water formation. These calculations instill confidence 
in the method employed here to detect areas favorable to bottom water 


formation. 
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7. Conclusions 

The stability formulas developed by Pollack are useful tools for 
the determination of vertical stability conditions. They are particular- 
ly adaptable to calculations using digital computers and offer results 
which are numerically consistent with results obtained with the less 
tractable formula of Hesselberg and Sverdrup. 

The calculation of density in situ and speed of sound in sea water 
using the formulas of Eckart and Kinsler-Frey extend the usefulness of 
the Pollack stability parameters in computer operation. These formulas 
offer accuracies comparable to those achieved through the standard ap- 
proach of using hydrographic tables. 

Based upon the stability calculations, bottom water was not formed, 
at least in the manner associated with the Nansen theory, in the vicin- 
ity of Station Mike during the period 1949-1953. Quanta of more dense 
and unstable water are shown by the calculations to exist in that area. 
However, these quanta do not penetrate the stable layer and reach the 
bottom. It is possible that they move laterally and downward, reaching 
the bottom waters in some other location. This cannot be investigated 
with the available data. 

The thickness of the bottom water layer at Station Mike undergoes 
aperiodic but definite variation. This is evidenced by the decreases 
in depth of the base of the stable layer. According to Crease (1965) 
the sill septh of the Greenland-Scotland Ridge is approximately 800 
meters. This is also the approximate normal depth of the base of the 
stable layer at Station Mike. However, during the period 1951-1953 
the base of the stable layer rose to a depth of about 400 meters on 


Zt 


six occasions. Cooper (1955) postulates that the current structure 
along the Greenland-Scotland ridge is such that it permits a pulsating 
rather than a continuous flow of water out of the Norwegian Sea, over 
the sill, and down into the Atlantic Basin. While this investigation 
does not relate to Cooper's theory it does indicate the possibility 


that pulses of rather large magnitude occur in the Norwegian Sea. 
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